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Abstract 

We study the baryogenesis via lepton number violation in the model of Anti-GUT. The origin 
of the baryogenesis is the existence of right-handed Majorana neutrinos which decay in a C, CP 
and lepton number violation way. The baryon number asymmetry is calculated in the extended 
Anti-GUT model which is only able to predict order of magnitude-wise. We predicted baryon 
number to entropy ratio, Yg = 1.46^'^ x 10 -11 , and this result agrees with experimental 
values very well. 

PACS numbers: 12.10.-g, 12.60.-i, 13.15.+g, 13.35.Hb, 13.60.Rj 

Keywords: Baryogenesis, Lepton number violation, Mass matrices 



January 2001 



* E-mail: hbech@alf.nbi.dk 
§ E-mail: yasutaka@nbi.dk 



1 Introduction 



The evidence of the neutrino masses and its mixing angles from the atmospheric and solar neutrino 
experiments Jl|, ||, ||], ||, |(| 0, [8|, ||] indicates a method to solve a challenging question in cosmology 
and particle physics: namely, baryon asymmetry of the Universe. This asymmetry cannot be 
explained with the "pure" Standard Model (SM) - with negligible B — L asymmetry - to the 
phenomenologically right magnitude of baryogenesis. In fact, the electroweak phase transition 
scenario for baryogenesis does not work very well in the frame of the SM with presently known 



lower bound of the SM Higgs mass [10]. Even for the minimal supersymmetric standard model 



(MSSM) this scenario is strongly disfavoured by baryon number wash-out at the electroweak phase 



transition ||ilj| . More detailed analyses of the MSSM [12] have been made and it is claimed that the 
MSSM is just consistent with baryogenesis in a very restricted region of parameter space requiring 
the right-handed stop to be lighter than the top quark and the left-handed stop heavier than 1 TeV. 

The model which we investigate in the present article - the extended Anti-GUT model - has as 
one of its characteristics that it coincides with the pure SM (at least) for energy scales below the 
see-saw neutrino scale, and so we have in this model no way to get the phenomenological baryon 
number unless we already have an B — L asymmetry prior to the weak epoch. 

However in the SM the neutrinos are massless due to the weak gauge symmetry or the conser- 
vation of the lepton number (L) and therefore the observed neutrino masses imply the existence of 
the new scale - the only well-known scales, i.e. weak, strong and Planck scales, alone are not able 
to provided this scale. A very suggestive mechanism is the the see-saw mechanism [[HI. 



In our extended Anti-GUT [14] model the see-saw mechanism is already built in as far as we 
assume the existence of the three right-handed neutrinos - Majorana neutrinos - in the range of a 
new scale in the order of 10 12 GeV, and the predictions of this model are very successful for the 
small mixing angle MSW E3| (SMA-MSW) scenario. 



Sakharov [16] has pointed out that a matter-anti- matter asymmetry can be dynamically gen- 
erated in an expanding Universe if the particle interactions and the cosmological evolution satisfy 
the three conditions: (1) baryon number violation, (2) C and CP violation, (3) departure from 
thermal equilibrium. The Anti-GUT model is a natural extension model of the SM in which C 
and CP are already violated, i.e. this model satisfies if combined with the standard cosmology all 
these Sakharov's conditions, therefore the baryon number asymmetry should be predictable with 
this model in the scenario that the three Majorana neutrinos are very heavy and that they violate 
the lepton number conservation during their out-of-equilibrium decays in the early stage of the 
Universe: Baryogenesis via lepton number violation [17]. 



The presently "observed" baryon asymmetry, the baryon density-to-entropy density ratio of the 
Universe, 

Yb = n B~ n B = (LQ _ 1Q) x 1Q -U (1) 

s 

is with the Universe evolution explained as consequence of the spectrum and interactions of the 
particles which break the linear combination of the baryon number (B) and lepton number, B — L, 
global symmetry at high temperature. 

In the following section, we should review briefly the Anti-GUT model to calculate the baryoge- 
nesis. Then, in the next section, we shall discuss the scenario of the baryogenesis via lepton number 
violation. In section 4 we present the results of our model. Section 5 contains our conclusion and 
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resume. 



2 The extended Anti-GUT model 

In this section we shall review briefly the extended Anti-GUT model [|l4j]. The extended Anti- 
GUT model is based on a large gauge group which is the Cartesian product of family specific gauge 
groups, namely, 

i= x lt2i3 (SMGiXU(l) B _ Lii ) , (2) 

where SMGi denotes 5{7(3)j x SU (2)j x U(l)i (SM gauge group), and i denotes the generation, i.e. 
each "proto-family" has a certain subgroup of the grand unification group, 50(10). This group in 
Eq. (H|) consist only of those representations that do not mix the different irreducible representation 
of the SM and is spontaneously breaking down to the SU(3) x 577(2) x {7(1) x U(1)b-l at the 
scale about 1 to 2 orders of magnitude under the Planck scale. The breaking is supposed to occur 



Table 1: All {7(1) quantum charges in extended Anti-GUT model. 
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by six Higgs fields which we have invented and denoted by the symbols S, W, T, £, x an d (f>B-L 
and their quantum numbers are given in Table 1. Finally the breaking of SU(2) x U(l) of the SM 
is broken by Weinberg-Salam Higgs field, (j) ws . (also its quantum numbers are found in Table 1.) 

We summarise here the vacuum expectation values (VEV) of the seven Higgs fields which the 
model contains: 



1) The smallest VEV Higgs field is the Standard Model Weinberg-Salam Higgs field, 0ws> with 
the VEV at the weak scale being 246 GeV/V2. 

2) The next smallest VEV Higgs field is also alone in its class and breaks the common B — L 
gauge group U(1)b-l, common to the all the families. This symmetry is supposed to be 
broken (Higgsed) at the see-saw scale as needed for fitting the over all neutrino oscillation 
scale. This VEV is of the order of 10 12 GeV and called 4> B -l- 

3) The next 4 Higgs fields are called £, T, W, and x an d have VEVs of the order of a factor 10 
to 50 under the Planck unit. That means that if intermediate propagators have scales given 
by the Planck scale, as we assume, they will give rise to suppression factors of the order 1/10 
each time they are needed to cause a transition. 

4) The last one, with VEV of the same order of the Planck scale, is the Higgs field S, which 
gives little suppression when it is applied, of the order of a factor l/\/2. 



The quantum numbers of the 45 well-known Weyl particles and additional three particles - 
Majorana neutrinos - are gotten from the requirement that all anomalies evolving U(1)b-l i, 
U(1)b-l,2, U(1)b~l,3 vanish strongly even without using Green-Schwarz anomaly cancellation 
mechanism jl8|| , i.e. the extended Anti-GUT model is an anomaly free model. These quantum 
numbers are also shown in Table 1. 

Now we can write down the mass matrices which are necessary to discuss the mechanism of 
baryogenesis in the following sections: the Dirac neutrino mass matrix and the Majorana neutrino 
mass matrix. These matrix elements were gotten using the technical corrections, factorial factor 
corrections, V # diagrams multiplying the mass matrix elements which take into account the possi- 
bilities of permuting the contributing Higgs fields. With this technical correction and the quantum 
charges of the Higgs fields the mass matrices, Dirac neutrino and the Majorana neutrino, are given 
by: 



D 



M, 




35SWT 2 £ 2 60^/14 S 3 WT 2 £ 3 
35S 2 WT 2 £ 2V3WT 2 
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2ws 3 x 2 e Vi5S x 2 ^ Vqs x C\ 
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SlEWT 2 
V6WT 



X 



(3) 



(4) 



We know neither the Yukawa couplings nor the precise masses of the fundamental fermions, but 
it is a basic assumption of the naturalness of our model that these couplings are of order unity and 
random complex in the Planck scale. In the numerical evaluation of the consequences of the model 
we explicitly take into account these uncertain factors of order unity by providing each matrix 
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element with an explicit random number \j - with a distribution so that its average (logAj,-) ~ 
and its spreading is 64%. Note that the random complex order of unity factors which are supposed 
to multiply all the mass matrix elements are not presented here. 



3 Fukugita and Yanagida scenario for the lepton number produc- 
tion 



The weak SU(2) instantons 1£] - sphaleron [20] - guaranteed the rapid exchanges of the baryon 



number and lepton number in which though B — L is conserved in the time of big bang, when the 
temperature was above the weak scale. But in our model we have the three right-handed neutrinos 
decaying in the L quantum number violating way, in fact also B — L violating way, at the time scale 
of the see-saw neutrinos. Therefore the baryon number violation condition of the first Sakharov's 
condition was effectively fulfilled. 

The assumption in our model that all the coupling constants and coefficients are of order of 
unity and random at the Planck scale, especially having random phases as far as allowed, implies 
not only C violation but also CP violation. 

Finally the third condition among the Sakharov conditions - out-of-equilibrium - comes about 
during the Hubble expansion due to the excess of the three type of the right-handed Majorana 
neutrinos caused by their masses. From these statements our model is seen to implement the 
scenario of Fukugita and Yanagida [17]. 



In the scenario favoured by our model the heaviest one among the three right-handed neutrinos 
turns out to gives the dominant contribution to the baryon or rather (B — L) quantum produced 
with the second becoming almost same order magnitude. Also it turns out that the average lifetime 
of this heaviest right-handed neutrino is of the same order as the Hubble expansion time so that 
we can count that a major part of these right-handed neutrinos first decay after inverse decays 
have essentially stopped. We shall justify and discuss these features of the scenario induced by our 
model in the next subsection. But first we shall describe the appearance of the baryon asymmetry 
taking for granted the mentioned assumptions so that we, (1) use only the heaviest right-handed 
neutrino and (2) assume it to live relatively long compared to the time needed before B — L is 
effectively conserved again. Really as we shall discuss more below, it is not the full B — L which 
is the important quantity but a special roughly speaking "third generation B — L" that is the 
sufficiently well conserved charge quantity. 



3.1 Conservation of the B — L-quantum charge 

A priori the excess of (B — L) quantum number risk to be diluted or washed out before the 
"accidental" (B — L) conservation of the SM sets in. It is therefore very important to argue for 
that such wash out does not take place. 

An important point is that it actually will turn out that if we only thought in terms of {B — L) 
one would estimate an appreciable wash out. However, we define a new "new charge" (B — L)^ 
which^ is washed away much more slowly. In the era until the lightest right-handed neutrino has 

1 Note that this (B — L)s is not exactly the same as the (B — L)z in our model, the latter a proto (B — L)s while 
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become so hard to produce that there are basically no more inverse decay processes producing it 
going on and also 2-by-2 scatterings are supposed negligible, we have effective conservation of the 
following charge (B — L)^: we define (B — L)$ to the baryon number minus the lepton number 
sitting on those leptons (or quarks) which are capable by collision with a Weinberg-Salam Higgs 
particle to produce, in resonance say, the heaviest of the right-handed neutrino, but not the two 
lighter ones. 

It the era when we can ignore diagrams involving the heaviest see-saw neutrino a "third gener- 
ation" quark or lepton, i.e. not coupling to the vertices N2 4> W s (■ or N\ 4> ws i cannot get converted 
in a (B — L)-violating way because the diagrams have to contain such vertices. 

The protection against dilution by N2 and N\ effects of the by N3 produced (B — L) hoped for 
is thus not relying on our model having a gauged (B — L)z but is a more general mechanism. 

The question of whether the {B — L) quantum number produced in excess by the decays of the 
heaviest Majorana neutrino will be preserved for the future is thus the question of whether the 
temperature falls so deep that this right-handed heaviest neutrino itself gets so hard to produce 
that we can ignore its inverse decay before all these heaviest right-handed neutrinos have decayed 
except for a fraction of order unity. The time needed to make the heaviest right-handed neutrino 
effectively unproduceable is the Hubble time corresponding to the temperature being equal to the 
mass of this see-saw neutrino. The crucial parameter to settle if this approximation of sufficiently 
slow decay is thus the ratio 



v _ M Planck ((M?)iM?)i 

Kj = — 



2H 



T = M > 1.66(^ s ) 2 8tt 5 1/2 Mi 



(i = 1,2,3) , (5) 



where Tj is the width of the flavour % Majorana neutrino, M ? - is its mass and 17*3 is the number of 
the degree of freedom at temperature Mj (see Eq. (7)). 



3.2 Baryogenesis and CP violation 



Now a right-handed neutrino, Nr, decay into a Weinberg-Salam Higgs particle and a left-handed 
lepton or into the CP conjugate channel. These two channels have different lepton numbers ±1. 
But, because of our random complex couplings, the partial widths do not have to be the same 
to next-to-leading-order perturbation theory. Defining the measure e ? - for the CP violation in the 
decay of the right-handed neutrino 

e * = r XT ' (6) 

1 N R .e + 1 NR .e 

where T N ^ e = Y<a,p T ( N Ri ^4>ws) and T Nr j = Y, a ^{N Ri -► l a 4p^ s ) are the N Ri decay 
rates (in the N Ri rest frame), summed over the neutral and charged leptons (and Weinberg-Salam 
Higgs fields) which appear as final states in the Nr. decays one sees that the excess of leptons over 
anti- leptons produced in the decay of one N Ri is just 

At high temperature (T>M^) equilibrium there were as many Majorana neutrinos per species 
as massless fermions, SM fermions. In the case (T^M^) all these Majorana neutrinos decay after 
the B — L violation is switched off. 



(B — L) s deviates by mixing angles from the first one. 
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To be able to calculate Baryogenesis we need to obtain the total number of effectively massless 
degree of freedom of the plasma, g*i, at temperature of the order of lightest right-handed neutrino 
(about 10 6 GeV), i.e., there are 14 bosons and 45 well-known Weyl fermions plus i Majorana 
particles: 



9* i 



T; 
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i = 2 
i = 3 



(7) 



here Tj denotes the effective temperature of any species j. When we have coupling as at the stage 
discussed between all the particles Tj = T. The entropy of Planck radiation with the degree of 
freedom j is 



2tt g*i 3 
45 



(8) 



Moreover, we should note here that due to the electroweak sphaleron effect, the baryon number 

ton number asymmetry Y 

Yb = a Yb-l = — —~r Yl 



asymmetry Yb is related to the lepton number asymmetry Yb by [21]: 

a 



with a 



a-1 
8N f + 4N H 
22N f + 13iVtf 



where Nf is the number of generations and Nh the number of Higgs doublets, this reads in the SM 
a = 28/79. 

Because Ki is not small, we have to expect a dilution effect for which we define the suppression 
factor Ki, i.e. we define Ki so that the resulting relative to entropy density, S{, baryon number 
density 



Y B 



y* % 



i=i 



A good approximation for Ki, the dilution factor, is inferred from Ref. [22, 23 1: 

0.3 



10<iQ<10 6 : 



1< ^ <10 : 



0< K <1 • 



Ki{hxKi 
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2Ki 
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6 



0) 

(10) 

(11) 
(12) 



Note that these dilution factors - we are taking it - contain the effect of the sphaleron processes 
we should use, instead of the Eq. (17) and Eq. (18), we took the following interpolating redefined 
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Figure 1: Tree level (a), self-energy (6) and vertex (c) diagrams contributing to heavy Majorana 
neutrino decays. 



dilution factor in the range 0^, Ki ^10: 

0< Ki < 10 : Ki = 1 (13) 

Since this dilution factor is smoother than the in the Ref. [^] defined one, and due to using 
order-of-one complex random factors in calculation of i"Q's, especially in our model, baryogenesis 
comes mainly from the third Majorana neutrino decay, i.e. the calculation of the K% have to be 
carefully performed because of if 3 ~ 1, therefore we should better use the newly defined smoothed 
out Ki in Eq. (19). 



3.3 CP violation in decays of the Majorana neutrinos 



The total decay rate at the tree level (Fig. 1 — (a)) is given by 



Nit 



+ r 



((M^)tM, 



4tt ' 



Mi 



(14) 



The CP violation in the Majorana neutrino decays, e^, arises when the effects of loop are taken 
into account, and at the one-loop, the only CP asymmetry of the vertex contribution comes from 
the diagram shown in Fig. 1 — (c). However there is an other contribution, from the wave function 
renormalisation (Fig. 1— (6)), which must be also taken into account and which gives typically same 



order amount as the vertex one 24, 25 



e,: 



^{^sYiiMPyMP), 



D ) 2 } 



(M}\ (Mf 



(15) 



where the function, f(x), comes from the one-loop vertex contribution and the other function, g(x), 
comes from the self-energy contribution. These functions can be calculated in perturbation theory 
only for differences between Majorana neutrino masses which are sufficiently large compare to its 
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decay widths, i.e. the mass splittings satisfy the condition, \M% — Mj\ S> |r» — Tj\: 

l + x~ 



/(*) 



1 - (1 + x)ln- 



x 



(16) 



g{x) = ^— . (17) 
1 — x 

4 Baryogenesis calculation in Anti-GUT model 

In this section we present the calculation of the baryogenesis at first numerically, and then in a 
crude way by seeking the dominating terms. 

The calculation goes by using formula fll5|), (16) and (17) taking as the mass matrix M®, 



the expression (3), with each matrix element further provided with an (independent of the other 
matrix elements) random complex number of order unity. Also the right-handed neutrino masses 



to be used as Mj (i = 1, 2, 3) in the functions / and g in formulas (|iq ) and (17) are (in principle) 
calculated by inserting random numbers of order unity into, in this case, the right-handed neutrino 
mass matrix @, in the way that each matrix element is again provided with an order of unity 
random (complex) factor. These random order unity coefficients are in the calculation taken as 
complex pseudo-random numbers. 

The quantity is then calculated 100,000 times and its logarithm, ln|ej|, is averaged over the 
different sets of random numbers. And the parameters which are the suppression factors S, W, T, £, \ 
and 4>b-l are taken for the best cases of our previous article fitting the neutrino masses and its 
mixing angles: 

(<j) ws ) = — GeV , {<j> B _ L ) = 2.74 x 10 12 GeV , (S) = 0.721 , 
(W) = 0.0945 , (T) = 0.0522 , (£) = 0.0331 , (x) = 0.0345 , (18) 

where the vacuum expectation values, except the Weinberg-Salam Higgs and {<P B -l)i are presented 
in the Planck unit. 

That is to say we really use 5, W, T, £, \ an d 4>b~l from one of the best fits including the 
factorial corrections to the charged lepton and quarks, while x an d 4>b-l are fitted to the neutrino 
oscillation data. Since our model unavoidable predicts the small mixing angle MSW solution, our 
calculation is not meaningful unless the SMA solution is right. But recently Barbieri and Strumia 



[27] have studied the neutrino oscillation fit using "non-standard analogies" method and have shown 
that the region of the SMA-MSW solution should be shifted in the direction of smaller mixing angle, 
so that it escapes the day-night effect exclusion by the Super-Kamiokande measurements. 



4.1 Results 



Our total baryogenesis is a sum of three contributions, one from each see-saw particle. The 
total baryon density-to-entropy density ratio is given by 



3 
i=l 



E6j 
y*i 



i=l 



(19) 



9 



where Yb % is the contribution due to the decay see-saw neutrino number i counted upwards in 
mass. 

Dominantly the signs of £3 and e 2 are strongly correlated but since the e3-contribution is diluted 
away it does not matter and the danger of cancellation is not there. Actually we included in the 
numerical calculation the correlation of signs effects correctly. In our case the contribution of Yb 3 
is dominant and the other contributions are negligible compared to Yb,3- This makes the scheme 
opposite to SUSY constraint ones. 

The numerical results of our best fitting case gives 

|e 3 | = 6.8 x 1(T 9 , K 3 = 1.06 (20) 
\e 2 \ = 6.0 x 10" 9 , K 2 = 4.29 (21) 
|ei| = 4.8 x 10~ 10 , K x = 19.8 . (22) 



Using the philosophy of letting order unity random numbers being given by a Gaussian distri- 
bution presented in the article p8| we estimate the uncertainty in the natural exponent for Yb to 
be 64 % • y/7 « 150 %. With these error estimations we get 

Y B = 1-46+^ x 10- 11 . (23) 
5 Conclusion and resume 

We have obtained with the cosmological fit, 

YB (experiments = (0-1 " 1.0) X lO" 10 , (24) 



Table 2: Typical fit including averaging over 0(1) factors. All quark masses are running masses at 
1 GeV except the top quark mass which is the pole mass. 





Fitted 


Experimental 


m u 


3.1 MeV 


4 MeV 


m d 


6.6 MeV 


9 MeV 


m e 


0.76 MeV 


0.5 MeV 


m c 


1.29 GeV 


1.4 GeV 


m s 


390 MeV 


200 MeV 


«v 


85 MeV 


105 MeV 


M t 


179 GeV 


180 GeV 


m b 


7.8 GeV 


6.3 GeV 


m T 


1.29 GeV 


1.78 GeV 


V us 


0.21 


0.22 


v cb 


0.023 


0.041 


V ub 


0.0050 


0.0035 


J CP 


1.04 x 10~ 5 


2-3.5 x 10" 5 



10 



excellently agreeing prediction of the baryon number-to-entropy ratio 



Y B = 1.46+^ x HP 11 , (25) 

from our extended Anti-GUT under use of parameters all already fit to either the charged lepton 
and quark spectra or to the neutrino oscillations. 

Even the "discrete fitting" of the precise choice of discrete quantum numbers of the seven Higgs 
fields (f>ws, S, W, T, £, % and 4>b-l of our model were fit already to the mass and mixing angle 
data. 

It should be remarked that our model predicts all fermion masses (neutrino mass square dif- 
ference ratio) and their mixing angles, including Jarlskog triangle, using above presented seven 
Higgs VEVs (The quarks and charged lepton mass spectra and their mixing angles are presented 
in Table 2.): 



Am% +30 o , . 

t ° =5.8 x 1(T 3 (26) 
A™ 2 tm -5 

-1-21 

tan 2 6» Q = 8.3 x 10" 4 (27) 
—6 

tan 2 (9 e3 = 4.3 +11 x 10" 4 (28) 

tan 2 atm = 0.97^ (29) 

i.e. our predictions are compatible with the MSW-SMA solar neutrino solution. Because of small- 
ness of the Cabibbo angle induces our solar mixing angle to be so small that it would stress the 
model drastically to seek to fit one of the series of large solar mixing angle fitting regions. 

All of the input parameters are, seven Higgs VEVs, already determined before the calculation 
of baryogenesis, so in this sense the baryon number-to-entropy ratio is pure prediction from our 
model! 

The number of measured quantities which are predictable with our model, quark and charged 
lepton masses and its mixing angles including the Jarlskog triangle area Jcp and also two mass 
square differences for the neutrino and the three of their mixing angles, is 19. Our model success- 
fully predicted all these quantities using only six parameters^]: the genuine number of predicted 
parameters is thus 13. But we have taken into the predictions the quantity tan 2 #13 for which 
CHOOZ has only an upper bound. 

It should, however, be admitted that our prediction is only order of magnitude-wise because we 
needed to use the assumption that all the coupling constants at the fundamental scale are of order 
unity. 



5.1 What do we learn? 



In the crudest approximation our agreement means that the general success of getting good 
baryon number prediction from the Fukugita-Yanagida scheme of having B asymmetry from the 

2 The VEV of Weinberg-Salam we have not counted as a parameter because of the relation to the Fermi constant. 
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see-saw scale works well with our specific model. This very good agreement of our prediction, 
of course, suggests that our specific model may carry some truth. It is first of all via the CP- 
violating parameter, e, that different detailed models can make their differences felt since even a 
supersymmetric model doubling of g* is only a factor 2 hardly distinguishable with our only order 
of unity accuracy. 

We expect the e which is an expression for the overall size of the Yukawa couplings - it is a loop 
effect relative to the tree diagram - for given mass splitting - for charged leptons and quarks to be 
sensitive to the number of effective charges in the model used for the mass protections. Indeed we 
expect more general suppression and therefore smaller e and thus baryon number if we imagined 
to have control over the ftj's - the bigger the number of charge types used. 

When our model then gives a good baryon number result we should expect that its number 
of charges species - effectively - used is roughly right. Since our model has in fact under some 
restrictions the maximal number of gauge charges the success of predicting e suggests that the right 
model should have a rather large number of charges species. It should, however, in this connection 
be reminded that since we have the field S giving only a tiny suppression the tree group we have 
"effectively" used is not the full Anti-GUT one but rather the subgroup of it obtained after breaking 
by S, SU(3) 2 x SU{2) 2 x U{lf. 

It should be also be mentioned that it was in our Anti-GUT put in, that the representations of 
the non-abelian subgroups SU(3) and SU(2) were given by a rule from those of the abelian, so that 
the non-abelian ones play no separate role. We did not even list the non-abelian representations in 
Table 1. 

Moreover, it should be stressed that our scenario is very different from the SUSY model one: 
in SUSY model or SUGRA there is the problem that gravitinos survive and cause an unacceptable 
mass density contribution unless their production does not occur due to late inflation p9[ . This 
makes it a problem to obtain the B — L from the see-saw neutrino decays and more problematic 
the heavier the right-handed neutrino used. This is why Buchmiiller and Pliimacher |^] have the 
lightest of the right-handed neutrino provide the baryogenesis. 

Since our model is a non-SUSY one we have no gravitino problem and thus we could equally 
well make use the heavier Majorana neutrinos as of only the light one. 

Acknowledgements 

We are grateful to W. Buchmiiller, CD. Froggatt, K. Kainulainen, M. Pliimacher and S. Lola for 
useful discussions. We would like to thank especially T. Yanagida for discussion of the parameter 
K. We wish to thank the Theory Division of CERN for the hospitality extended to us during our 
visits where part of this work was done. H.B.N, thanks the EU commission for grants SCI-0430-C 
(TSTS), CHRX-CT-94-0621, INTAS-RFBR-95-0567 and INTAS 93-3316(ext). Y.T. thanks the 
Scandinavia-Japan Sasakawa foundation for grants No. 00-22. 
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